Effects of Density of States Asymmetry on Boundary Resistance of 
Ferromagnetic-Nonferromagnetic Metal Interface 
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Asymmetry of conduction of electrons in the spin majority and minority bands of a ferromagnetic 
metal is well known to produce spin accumulation at the boundary with a normal metal when a 
current is injected through the interface. However, little emphasis has been put on the density of 
states (DOS) asymmetry in diffusive ferromagnetic-nonferromagnetic multilayers. We found that if 
DOS and conduction asymmetry differ, the electric potential in a ferromagnet falls off to its bulk 
value exponentially on a scale of the spin diffusion length. Therefore the boundary contribution 
to resistance dramatically depends on whether the potential difference is measured close to the 
interface or farther than the spin diffusion length from it. This result is not altered by taking 
surface resistance or spin-flips on the intermetallic boundary into account. Explicit answers for 
common multilayered spin-valve structures are given. 

PACS number(s): 72.25.Pn, 72.25.Ba, 73.40.Cg, 73.40. Jn 



I. INTRODUCTION 

The possibility to exploit electron spin to manipulate 
currents in electronic microdevices has generated wide 
interest in the rapidly expanding field of spin electronics 
01 QUI- Technically, the itinerant electron spin may be 
put to use if the spin-flip scattering, intermixing the spin- 
up and spin-down channels, occurs at a rate r^ 1 which is 
slow compared to the inverse mean free time r _1 which 
characterizes the dissipative processes in charge of electri- 
cal resistivity. Separation of electrons in two channels is 
a key element of the physics of a central spintronics phe- 
nomenon of giant magneto-resistance (GMR) of metallic 
ferromagnetic-nonferromagnetic metal multilayers. 

The two-channel description of diffusive electrical 
transport from a ferromagnetic (FM) to a nonferromag- 
netic metal was given in Q and on a thermodynamic 
basis also in Q. Due to the exchange splitting of the 
Fermi surface in a FM metal the Fermi energy corre- 
sponds to different parts of the band structure for spin- 
up and spin-down electrons, and therefore the two types 
of electrons have different densities of state iV-^ and con- 
ductions o-fjj,. Because of the conduction asymmetry 



cr T + a i 



(1) 



an electrical current passing from a FM to a nonferro- 
magnetic metal is polarized and injects magnetization. 
This effect shows up as an additional contact resistance. 

The size l s d of the region of the spin accumulation is 
set by the equilibrium between spin injection through the 
interface and the spin-flipping processes in the bulk of the 
paramagnet. The spin diffusion length l s d relates to the 
mean free path A as Z s d/A ~ (t^i/t) 1 / 2 . An analysis of 



the two-channel model was given by Valet and Fert @ 
on the basis of the Boltzmann transport equation in the 
limit l s d 3> A. They have also applied the model to the 
geometry of the GMR with current perpendicular to the 
plane (CPP-GMR). A significant advantage of their ap- 
proach is the ability to describe measurable quantities in 
detail in terms of parameters of the theory, which allows 
to extract the latter from experiment [a,|9, 10]. 

The formulas in [fj, however, have been derived in a 
flawed assumption that while the conductivities a^,[ of 
the majority and minority spins are different, the corre- 
sponding densities of states (DOS), N^i, are equal. But 
the difference in N^ i is essential as it leads to a jump of 
the electric potential on an intermetallic interface, which 
also was not observed in the later work jj. Neglect of 
the difference in results in an incorrect variation of 
the electric potential with distance from an intermetallic 
interface, and changes, inter alia, the boundary contri- 
bution to resistance. 

The aim of the present paper is to provide expressions 
for boundary resistances of common multilayered systems 
corrected to allow the measurements of potential drop to 
be performed at any distance from interfaces, not neces- 
sarily ^> l s d- Correlative results for an array of domain 
walls have been obtained in Hfl. In Sec. HD a review 
of the two-channel diffusive transport theory is given, 
whereas in Sec. IHII we demonstrate our results for a single 
ferromagnetic-nonferromagnetic interface and for typical 
spin-valve structures. Sec. IIVI summarizes the findings 
and provides a brief discussion. Finally, an appendix with 
details of our approach to calculation is presented at the 
end of the paper. 



II. BASIC EQUATIONS 
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Starting from the Boltzmann equation, Valet and Fert 
0| showed that in the limit i s d ^ A one can describe 
transport in the system introducing electrochemical po- 
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tentials p p for the two spin channels p =T> j as sums of 
chemical potentials (accounting for the kinetic energy of 
electrons) and the potential energy. For small deviations 
from equilibrium the chemical potential is given by the 
non-equilibrium electron density n p divided by the den- 
sity of states N p on the Fermi surface: 



p p = rip/Np - eU. 



(2) 



Here — e is the electron charge, U is the electric poten- 
tial, and the condition n-f = — is fulfilled owing to 
electroneutrality. 

We now briefly review the equations of two-channel dif- 
fusive transport to establish notations. Electron current 
j p in the channel p is a sum of a diffusion component and 
a contribution from the electron drift under the influence 
of an electric field E = —dU : 



3p = -{-e)D p dn f 



crpE 



a-pdpp/e, 



(3) 



where we used the Einstein relation a p = e 2 N p D p be- 
tween the conductivity a p and the diffusion constant D p . 
In the bulk the excess electron density n-f is zero, the 
electrochemical potentials for the two channels are equal: 
/if = /if, and the current polarization is a JIJ. Relax- 
ation of non-equilibrium polarization of current near a 
boundary between two materials due to spin-flip pro- 
cesses is described by the continuity equations 



dip = en p /r n . 



(4) 



It is convenient to introduce the divergence p s = /if — 
Pi between the two electrochemical potentials. It is non- 
zero only in the presence of spin accumulation: from <|2j 
one finds that it is related to rif via 



p s = n 1 (N- 1 + N' 1 ). 



(5) 



Substituting (JjJ into <EJ and J3J we find that p s satisfies 
the diffusion equation 



d 2 p s = p s /l 2 d , 
where the spin diffusion length 



'sd 



n; 



(6) 



(7) 



Equation 101 by means of relation J&J describes the expo- 
nential decay of non-equilibrium spin density away from 
a boundary. 

In order to find currents j p in each channel and the 
electric potential U we introduce 



Mo 



erf/if + &iPi 



(8) 



where a = <7f + er^ is the net conductivity. Differentiating 
/to yields j = crdpo / e , where j = jf + Jx is the net current 
being conserved throughout the sample (dj = 0). The 
spatial distribution of po is found from d po — 0, hence 
in one dimension 



Po(x) — const +ejx/a. 



(9) 



The spin current j s = jf — j j may then be found from the 
known distribution of p s : 



3s = aj + f(l - a 2 )dp s /2e. 



(10) 



From the definition <|2J we find that the spatial distribu- 
tion of U is given by 



-eU = 



+ N l p l 
N 1 +N l 



Mo 



(11) 



where we have introduced the asymmetry of the density 
of states 



(3- 



N t + Ni 



(12) 



Equations (HJ, ©, flOll - and ltTT|l completely describe 
electron transport inside a layer. 

At distances 3> l a d from boundaries, when the diver- 
gence p s decays to zero, the two electrochemical poten- 
tials p^.i = po + |(il — a )Hs as well as the electric 
potential blU collapse to pq. I.e. po is the value of 
the electrochemical potential that would be realized were 
there not non-equilibrium distribution of spin. However, 
if measurements are done closer than l s d to an interface, 
one has to take into account that the non-equilibrium 
divergence of /if 4 has not yet subsided. 

To find charge and current distribution through the 
whole multilayer sample, Eqs. ©, illOll . and 11 111 
have to be supplemented by boundary conditions. In the 
absence of an interface resistance and spin-flip scattering, 
the appropriate boundary conditions are the continuity 
of pi 4 and jfj 0| . The continuity of /if 4 shows that po 
and U are discontinuous at intermetallic interfaces with 
the jumps 



No = \ [a]p s 



-e[U] = W]p 



(13) 



where [a] is the jump in the conduction asymmetry (QJ 
and [P] is the jump in the asymmetry of the density of 
states. The discontinuity in po appears in transport mea- 
surements farther than i s d from an interface as a bound- 
ary contribution to resistances 



Rio 



No/ej = [a]p s /2ej. 



(14) 



Discontinuities of the electrochemical potentials and 
the spin current may be introduced to account for in- 
terface resistance and spin-flip scattering. One then has 
to allow for the jumps of /if 4 in l|T3jl . which affects the 
answer H14I1 . But the inclusion of new parameters for 
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A Single Ferromagnetic-Nonferromagnetic 
Metal Interface 
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Figure 1: Sketch of spatial variation of potentials at the 
boundary between ferromagnetic and normal metals. 



interface resistance and spin-flip scattering renders the 
formulas very cumbersome, not changing the principal 
point of the paper. Therefore, we have delegated the dis- 
cussion of the interface resistance and spin-flip scattering 
to Appendix. 



III. RESULTS 

Equation H ill presents the spatial variation of the elec- 
tric potential, provided /io and /i s for the problem in 
question are found. As revealed by its second term, if 
(3 — a ^ 0, spin accumulation near an intermetallic in- 
terface causes U to collapse to its bulk asymptotic value 
—/io/e on a spin diffusion length scale. Hence, the form 
of the boundary contribution to resistance has to be mod- 
ified dependent on how close to the boundary the mea- 
surements take place. 

Anyhow, since the jump — e[U] at an interface as given 
by the second of Eqs. l(T3"j) differs, generally speaking, 
from the jump of [/x]o, if measured in the immediate vicin- 
ity to the interface, instead of 11411 the interface resistance 
would rather be 



Rio = ~[U]/j 



'K/2ej. 



(15) 



Note that the jumps — e[U] and [/i]o not only differ in 
absolute value, but may in principle have opposite signs. 

At intermediate distances from an interface the bound- 
ary resistance must be obtained from ltTT|l : R X1X2 — 
(U(x\) — U(x2)) /j by a complete solution of the problem 
for the whole multilayer structure. 



The solution is straightforward in the simplest case of 
an isolated boundary at the origin between two differ- 
ent materials with constants denoted by the two sub- 
scripts 1 and 2. We here have decided to illustrate our 
claim with more general than simply a ferromagnetic- 
nonferromagnetic interface example at an expense of only 
slightly more complex expressions. 

Equating spin currents ltTn|l we find [i s at the origin and 
reestablish the Valet-Fert expression for the boundary 
resistance: 



Rloo — 



(ai 



a 2 ) 



Ai + A2 



(16) 



In this expression and below we use the shorthand 
notation 



A = <r(l - a 2 )/l sd . 



(17) 



For arbitrary x\ < 0, x 2 > the resistance between the 
two points is 



Rxia 



\xi\/ai + x 2 /a 2 + -R/12, 



(18) 



where the first two terms are the series resistances of the 
two materials over the lengths x\ and x 2 respectively, 
and the interface contribution is 



-R/12 = Rloo 1 



(ai - (3!)e xi ^ dl - (a 2 - (3 2 )e~ X2 / l - d2 



ai — a 2 

This expression exponentially decays from the value 

lh-02 



Rio = Ric 



Q'l — Q'2 



(19) 
(20) 



very close to the boundary to i?/oo farther than l s d for it. 

Note that if one of the materials is a nonferromagnetic 
metal, then the corresponding exponential term in l|19ll . 
which describes decay of the electric potential in it, is 
zero. This is because in a normal metal the electric po- 
tential 111 111 falls onto the asymptote right after the 
interface, not exponentially with the distance from it. 

Position dependence of the electric potential ~eU and 
of the electrochemical potentials /io, /if, /i|, Ms f° r an 
isolated boundary between a ferromagnetic (ap = 0.5, 
(3-p = —0.5) and a normal metal (q;n = An = 0) is plotted 
in Fig. Note that because @f is taken negative, the 
jumps — e[U] and [/x]o have opposite signs. 

Allowing for interface resistance and spin-flips compli- 
cates the formulas 11 fill . Ill9|l . and l!20t without changing 
the qualitative conclusion of the difference between Rj 
and Rioo- Expressions for boundary resistances in this 
case are provided in Appendix. 
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Figure 2: Sketch of a spin-valve multilayer structure, with, 
supplying wires (two bottom structures) and without them. 
In each case parallel and antiparallel configurations of the two 
ferromagnetic layers are depicted. 



where 0n = tN/4dN- In the limit of large £n 3> 'sdN 
this expression turns into iflfi)) for an isolated boundary. 
When the thickness of the paramagnetic layer % — > 0, 
we obtain that Rioo vanishes as Rioo ~ a 2 t^/2a^i. 

In what follows we omit the series contributions to re- 
sistance like the first two term in l|18ll . For any structure 
in Fig. [2] these contributions are trivial. We thus sim- 
ply cite the boundary part of the total resistance of a 
symmetric trilayer between two probe points xl < and 
Xr > t N : 

Rirl = Rioo \ 2 - (e XL / l ° d * + e (*N-* R )/*sdF 

(22) 

For the measurement points very close to the boundaries 
xl — and xr •, = £n reduces to 



Rio = 2(3R Ioc /a. 



(23) 



For an antisymmetric trilayer without the supplying 
electrodes (the second structure from the top in Fig. |2j, 
when the right ferromagnetic layer has the opposite con- 
duction asymmetry —a, the Valet-Fert boundary resis- 
tances are 



B. Spin-valve Structures 

In the end, we cite formulas for a simplest spin- 
electronic device — a trilayer 'spin- valve': a (thin) layer 
of paramagnetic metal of thickness £n sandwiched be- 
tween two ferromagnetic layers with parallel or antipar- 
allel orientations of the magnetic momenta. We consider 
both the cases with and without supplying normal metal 
electrodes (see Fig. |2J. In the first case we assume the 
two ferromagnetic layers to have equal thickness tp. 

In case of a symmetric trilayer without the supplying 
electrodes (the uppermost structure in Fig. IKil the Valet- 
Fert boundary resistances of the two interfaces are 



Rio 



or 



Af + An coth #n /2 ' 



(21) 



Rioo = 



A F + A N th N /2' 



(24) 



When £n - * the resistance tends to a finite value 
Rioo = a 2 Z s dF/(l — ck 2 )ctf, thus justifying the term 'spin- 
valve'. The boundary contribution to resistance between 
two arbitrary points xi < and xr > £2 is 



Rirl = Ric 



and tends to 



(t-^—XR_)/l BC \F _j_ 



■0 



a 



Rio = 2(a - P)Rioo/a 



(25) 



(26) 



when xl — and xr ; = £n- 



For a symmetric five-layer with two supplying normal-metal electrodes (the third structure from the top in Fig. |2J 
the boundary resistances of the two external and the two internal interfaces are respectively 



Rle 
Rli 



An coth N /2 + A F th F /2 



A| + A N A f coth F (1 + coth N /2) + coth N /2 ' 
, A N + A F th0 F /2 

A| + A N A f coth F ( 1 + coth N /2) + A^ coth N /2 ' 



(27) 
(28) 



For an antisymmetric five-layer (the lowermost structure in Fig. |2j 



2 A N th0 N /2 + A F th F /2 



Rleo ° a A| +A N A F coth0F(l + th0N/2) + A2 J th0N/2' (29) 

d 2 A N + A F th F /2 

A| + A N AFCoth0F(l + th0N/2) + A 2 J th0 N /2' 1 ' 
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The boundary contribution to resistance between two arbitrary points xl < and xr > % + 2ip is for both 
symmetric and antisymmetric five-layered structures simply the total boundary resistance of all the interfaces Rirl = 
2(Riioo + Rieoo) since in the normal-metal supplying electrodes the electric potential coincides with its bulk asymptote 
— /ito/e everywhere after the interface. 



IV. CONCLUSIONS 

In this brief study we focused on the effects of the 
density of states asymmetry in the diffusive transport 
through multilayered ferromagnetic-nonferromagnetic 
metal structures. 

Whereas conduction asymmetry brings about spin ac- 
cumulation in a layer of thickness of the order of the 
spin diffusion length, and determines the overall drop of 
the electric potential in this layer (apart from the triv- 
ial series resistance contributions), the DOS asymmetry 
results in a jump of the electric potential directly at an 
intermetallic interface. In a standard Valet-Fert theory 
only the boundary contributions to resistance Rioo from 
the overall potential drop in the spin accumulation layer 
are calculated. These boundary resistances emerge only 
when measurements are done farther than the spin dif- 
fusion length from interfaces. 

However, the series resistances that build up over the 
distance ~ Z s( j between the measurement points turn out 
to be larger than the boundary resistance Rioo ■ To iden- 
tify the part of the resistance from the current conver- 
sion process on the boundary it is advantageous to attach 
probes closer to the interface. But in that case it should 
be remembered that because of the asymmetry of the 
density of states the interface contribution to the resis- 
tance differs from i?/oo- 
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conditions of continuity of the electrochemical potentials 
and of the currents apply. The form lf33~fl for a solution 
in each layer with respective Z s d relies on the values of 
/i s on the boundaries, which thus need to be put equal 
on the two sides of each interface to satisfy the bound- 
ary condition of the continuity of /i s = /if — /ij. The 
second boundary condition for the electrochemical po- 
tentials leads to the expression l|13ll for the jump of /io, 
which together with solution in each layer with re- 
spective a yields answer for the spatial distribution © 
of fi . 

Substituting into 111 PL we obtain spin currents at 
the boundaries: 



where 



Js[a) 



A 



aj + (-Bfi s (a) + Afi s (b))/2e, 

aj + {-Afi s (a)+Bn s {b))/2e, (32) 



sinh ■ 



B = A coth ■ 



(33) 



In case of a semi-infinite interval, when either b — > oo or 
a — ► — oo, fi s on the corresponding end tends to zero and 
A = 0, B = A. 

Equating 132t on both sides of each interface, we ob- 
tain an algebraic system of equations on the interface 
values of /i s . Once found, these values then give bound- 
ary resistances iTHI) and flU - In this way we derived Eqs. 
CSJ, (23, (23, (21, and (23, (Si- 
Interface resistance as well as spin-flip scattering at 
an interface may be taken into account by appropri- 
ate boundary conditions Q. Non-zero surface resistance 
leads to discontinuity of electrochemical potential pro- 
portional to the current in the respective channel and 
spin-flip scattering induces discontinuity of the spin cur- 
rent proportional to the average value of /i s : 



GENERAL SOLUTION OF THE EQUATIONS 

In this appendix we expound the general approach we 
used to calculate the spatial distribution of potentials in 
a multilayered structure. The diffusion equation © in 
an interval a < x < b in one dimension has a solution 



/i s (a) sinh ^-j 2 



-A*.(6) sinh^ 



sinh 



(31) 



that falls off exponentially on a length scale i B d from the 
values /i s (a) and /t s (&) at the boundaries a and b. 

First consider the boundaries with no interface resis- 
tance or spin-flip scattering. In this case the boundary 



Np = er p {j} 



pi 



[?']. - Ms, (34) 

er sf 



where {. . .} denotes an average of a discontinuous quan- 
tity between the values on the two sides of a bound- 
ary. The surface resistances r p in each channel may be 
parametrized by the total resistance 4r and the surface 
resistance asymmetry 
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(35) 



Solution 1-VH of Eq. inside a layer then obviously 
still holds, if /i s (a) and /i s (&) are understood as values 
of /i s on the internal side of each boundary. Substitut- 
ing this into expression lllOt yields the values of the spin 
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current on the internal side of each boundary, which may 
be combined with the boundary conditions (I34T) to give 
an algebraic equation on the interface values of ^ s . How- 
ever, it is easier to rewrite Eqs. and II32H in terms 
of the average values {/i} s and {j} s to produce: 

{j} s (a) = a r j + {-B r {n} s (a) + A r {LL} s {b))/2e, 
{j} s (b) = a r j + (-A r {n} s (a) + B r {n} s (b))/2e (36) 



in analogy with p2|| . Here parameters of interface re- 
sistance r, 7 and of interface spin-flip scattering r s f now 
enter only through the renormalized constants a r , A r , 
and B r : 



A.f 
B r 



a - (A + B)jr/2 
l + (A + B)r/2 ' 
A(l-r/2r s{ ) 
(l + Br/2) 2 - (Ar/2) 2 ' 
B + (B 2 -A 2 )r/2 + {l- 



Br/2)/r s 



(l + Br/2) 2 - {Ar/2f 



(37) 

(38) 

(39) 
I 



The formal analogy between l-'Uit and 13211 allows to 
reuse the system of algebraic equations on the interface 
values of fi s with the substitution of {^} s as variables 
and Q37J1 - H39H as parameters. However, the expressions 
for the jumps of /io and U at the boundaries are more 
complex than l(T3)l : 

Mo - er((l + {a} 7 )j + (W + 7){j} s ) + |W{^ft) 
-e[U] = e r((l + {/3} 7 )j + ({ ( S} + 7 ){i} s ) + I[/3]{^J.) 

Hence the boundary resistances 111411 . Hot will also be 
different. 

As an example and to give the reader a feel of the com- 
plexity of the formulas that arise when surface resistance 
and spin-flips are taken into account we cite our answers 
for a problem considered earlier in the text with simple 
boundary conditions of continuity, i.e. a single boundary 
at the origin between two materials. 



The boundary resistances farther than Z s d from the interface and close to it are respectively 



Rlao 

Rio 



r(l - 7 2 ) + 
r(l -7 2 ) + 



[a] 2 + r (Ai(a 2 + 7? + A 2 (ai + 7?) + 2r sf 1 ( 7 + {a}) 2 
Ai + A 2 + A 1 A 2 r + r sf 1 (2+ (Ai + A 2 )r/2) 
, , , [a}[/3] + r (Ai(a 2 + 7 )(/?2 + 7) + A 2 («i + 7) (A + 7)) + ar^ 1 ^ + {a})( 7 + {/?}) 



Ai + A 2 + AiA 2 r + ^(2 + (Ai + A 2 )r/2) 
The dependence of the boundary resistance on the distance from the interface is given by 



(42) 
(43) 



where 



2ej 2ej 

a] + r (Ai(a 2 + 7) + ({"} + 7)/4r sf ) 



Ms(+0) = 2ej 
H.(-Q) = 2ej 



Ai + A 2 + AiA 2 r + r sf 1 (2 + (Ai + A 2 )r/2) ' 
[a] - r (A 2 (ai + 7) + ({a} + 7 )/4r sf ) 



Ai + A 2 + AiA 2 r + r^(2 + (Ai + A 2 )r/2) 
are the values of /i s on the two sides of the boundary. 



(44) 

(45) 
(46) 
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